Rotational motion of a droplet induced by interfacial tension 
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Abstract 

Spontaneous rotation of a droplet induced by the Marangoni flow is analyzed in a two-dimensional 
system. The droplet with the small particle which supplies a surfactant at the interface is consid- 
ered. We calculated flow field around the droplet using Stokes equation and found that advective 
nonlinearity breaks symmetry for rotation. Theoretical calculation indicates that the droplet spon- 
taneously rotates when the radius of the droplet is an appropriate size. The theoretical results were 
validated through comparison with the experiments. 
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I. INTRODUCTION 



Recently, self-propelled motion of both biological objects and nonbiological objects are 
extensively studied to find general aspects in self-propelled systems [1]. Collective behaviors 
such as nonequilibrium order-disorder transition of direction of motion are candidates for 
such general characteristics [2]. One of the best experimental systems to investigate gen- 
eral aspects is self-propelled objects driven by interfacial tension [3, 4]. Actually, several 
results on collective behavior of the objects driven by interfacial tension have been reported 
recently [5]. This is because this kind of systems usually consists of a small number of ele- 
ments and allows us to perform experiments under broad parameter space. Despite of this 
advantages, however, mechanical analysis for above droplet motion is still lacking. 

The above-described droplet motion driven by interfacial tension is observed even under 
isotropic conditions. Although a droplet itself has no asymmetry, symmetry of interfacial 
tension is broken through the advective nonlinearity. Actually, Toyota et al. reported 
spontaneous translational motion of a circular droplet [4]. Recently, theoretical analysis 
for such translational motion has been reported and the nonlinear effect which breaks the 
symmetry has been made clear [6]. Not only spontaneous translational motion but also 
spontaneous rotation has also been demonstrated by Takabatake et al. experimentally [7]. 
In this system, a droplet rotates even though the system is mirror symmetric. Although this 
system was analyzed using Langevin equation in the previous study, mechanical analysis 
is lacking and the physical mechanism of symmetry breaking is not clear. In this article, 
we analyzed the spontaneous rotation of a droplet using an advection-diffusion equation 
coupled with Stokes equation. We found that the droplet can rotate when the droplet has 
an appropriate radius. To validate our model, theoretical results were compared with the 
experimental ones. 

II. MODEL EQUATION 

We consider a circular undeformable oil droplet with a radius, R, in an aqueous phase 
in a two-dimensional space, as shown in Fig. 1(a). A small solid particle is attached at the 
interface and a surfactant that reduces interfacial tension is supplied to both the inside and 
the interface of the droplet from the particle. In the droplet, the surfactant is decomposed. 
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FIG. 1. (Color Online) (a) Schematic diagram of the model in this article, (b) The flow field 
when the interface concentration field is c{6). The images are stream line around the droplet in 
the frame where the particle and the center of mass of the droplet are fixed. The red line is the 
oil-water interface and the white line is the particle, (b-1) is the flow field when c(9) = 1 + 0.1 cos 9 
(symmetric concentration field) and (b-2) is the flow field when c{6) = 1 + 0.2 cos 9 + 0.1 sin 9 
(asymmetric concentration field). 

As a result of these process, there appears inhomogeneity of interfacial tention, and the 
droplet is driven by the Marangoni flow induced by the interfacial tension gradient. In a 
laboratory frame, the particle can freely move along the interface. Without loss of generality, 
we may also consider a frame in which the particle and the center of mass of the droplet are 
fixed. Polar coordinates are introduced with the center of mass of the droplet regarded as 
the origin. The position of the particle is fixed at r — R and 6 = 0. With the assumption 
of low Reynolds number, the velocity field is calculated using Stokes equation as follows, 



- Vp + r]V 2 v = 0, 



(1) 



where v is the velocity, p is the pressure field and r\ is the viscosity coefficient. In our model, 
v satisfies the incompressive condition, V • v = 0. Due to the interfacial gradient, there is 
the following stress jump at the interface [8]: 



r=R 



(o) 
r0 



r=R 



1 <9 7 

Rd6' 



(2) 



where cr is the stress tensor, the superscripts "(i)" and "(o)" correspond to the in- 
side and outside of the droplet, and ^(8) is the interfacial tension between two flu- 
ids. Here, 7 linearly depends on the interface concentration of the surfactant, c{9) = 
c° + J2m=i {c c m cosm6 + c^sinm^), as 7(6?) = 70 - kc(9). 

Assuming that the surfactant supplied from the particle is distributed only at the interface 



due to fast decomposition in the bulk, the time evolution of c{6) can be described as 

dc u dc D d 2 c ..... , oX 

where D is the diffusion constant at the interface, a is the decomposition rate, and na/3 is 
the rate of supply. In Eq. (3), the stationary flow field satisfying Eq. (1) is used as u, which 
is determined by c(6). 



ILL SOLUTION OF THE TWO DIMENSIONAL STOKES EQUATION 

In this section, v satisfying Eq. (1) is calculated under given concentration field at the 
interface, c{6). In the polar coordinate, v is described as 

v = v r (r, 6)e r + v e (r, 9)e e , (4) 

where e r and e 9 are the unit vectors in the radial and angular directions, respectively. Due 
to incompressibility, 

V-v = -—(rv r ) + -— = 0; 5 
r or r do 

therefore, using the stream function, the solution of Eq. (1) can be described as 

Taking rotation of both sides of Eq. (1) and substituting Eqs. (6) and (7), we can obtain 

V 2 V 2 ^ = 0, (8) 



where 



V= = i— ( r —) + i *L. (9) 

r dr \ dr I r 2 d8 2 
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The solutions inside the droplet, and the outside the droplet, are 




The boundary conditions of flow field at the interface are v. 



(o) 



(o) 



(i) 

r=R 9 

free condition, 



r=R 



(») 



and 



r=_R 



r=R 



since the droplet shape is fixed. Since there is no external force, a torque 



= Jdl (n.^| r Jxr = 0, 



(12) 



and a force free condition, 

f = J dln.<r^\ r=R = 0, (13) 



have to be satisfied, where I is the distance along the interface. Using the conditions de- 
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scribed above and Eq. (2), the coefficients are calculated as 



n(i) _d(o) 
n ~ n 5 



C =0, 

B{ =B\ = 0, 

A\ = - Cf } = 
A\ = - Cf = 



(m^0,l), 
(m^0,l), 



KC S m R 



(14) 



where k = k/ {2 (77^ + 77^)}. From Eqs. (6), (7) and (14), flow profiles are obtained as 



,(°) 



,0) 



,(0) 



,0) 



E 

m=l 



E 



cos m# + sin m#) , 

m+l / \ m— 1 



(*)"" " (*)' 



m=l 



(c^ cos m# + c s m sin m#) , 



2-^-r 



i? 2 



-E 



«[(m-2) _,„(£) 



m+l 



m=l 



(c^ sin m6 — c s m cos m#) , 



2^r 



R 2 



-E 



K [-(m + 2) + m (J) 



r \ m— 1 



m=l 



(c^ sin m# — c s m cos m#) , 



(15) 
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and the tangential flow at the interface (r — R), u(9), is obtained as 

u{6) = + « X (c c m sin m# - cos m9) . (16) 

m=l 

As the particle is a solid, the flow velocity on the particle is (no-slip condition). As- 
suming that the particle size is small enough, this condition can be described as u(0) = 0. 
Using this condition, is determined as 



kR 2 c s 



= ~^rJ2f- ( 17 ) 

m=l 

At infinity, = —k (c\ cos 9 + c\ sin6>) /2 and vj^ = — — k (— c\ sm.9 + c\ cos 9) /2, 
which means the droplet translational velocity (V) and the angular velocity in the laboratory 
frame (Q) are 

V =K^e x + K^e y , (18) 



,s 



m=l 

e x and e y in Eq. (18) are the unit vectors in the direction of x axis (9 = 0) and y axis 
{9 = 7r/2), respectively. The calculated flow fields are illustrated in the case of translational 
motion and in the case of rotational motion using Line Integral Convolution [9] in Fig. 1 (b). 
It is noted that the rotational speed of a droplet is proportional to the sum of the Fourier co- 
efficients corresponding to an antisymmetric concentration field about the anterior-posterior 
axis. 

IV. ANALYSIS OF THE MODEL 

Using Eq. (3), we calculate the velocity and the angular velocity of the droplet in this 
section. Since higher Fourier modes of c decay quickly due to diffusion, we neglect c c n , and 
c s n for |n| > 2. Letting Cq, c\, and c\ be replaced by X, Y and Z, Eq. (3) is described as 

^ = - a(X - Ac) + 1 {F 2 - (y ) 2 + Z 2 } , (20) 

%=X,<Y-Yo)-±*, (2D 
^ =^Z + ^(Y- Y ) Z. (22) 



7 



Here, X = Kf3 2 p 2 /R (p 2 + 1) +/3/2, and Y = /3p 2 / (p 2 + 1) are the values of X and Y at the 
fixed point at which Z = 0. Here, p = R/1& is the radius normalized by the characteristic 



respectively, where L = / a // d is the ratio between another characteristic length, / a = 
K/3/a, and Id- l& is the length of advective transportation during 1/a since u ~ kc\ ~ 
K/9p 2 / (p 2 + 1) ~ k[5. If \ z is positive, the fixed point, (X , Y , 0), is unstable, which means 
a finite angular velocity of a droplet obtained from Eq. (19). 

The phase diagram of the motion can be illustrated as shown in Fig. 2 (a) using two kinds 
of motion: translational motion and rotational motion. To realize rotation, it is necessary 
that the advection term in Eq. (3) is dominant; that is, L ^ 1 is needed. In fact, \z can 
be positive only when L > 16-^/3/9 ~ 3.08. When -R ~ (p ~ 1) 5 the concentration 
field is almost uniform due to diffusion, so that flow is too weak to generate rotation. On 
the other hand, when R ^ Z a (p <; L ), the uniform decomposition term in (3), — ac, 
becomes dominant and the concentration field is almost uniform. Therefore, only when 
L is large enough and R is between R^ ps Z d and R^ w Z a , can the droplet spontaneously 
rotate due to the coupling between the mirror symmetric concentration field (F ) an d mirror 
antisymmetric concentration field (Z) through the rotational flow field. Since there is no 
other bifurcation, there are only a pair of stable fixed points with non-zero Z, which are 
calculated as X = 3a/ a /2ft, Y = —R\ y /k and Z = ±R\ / —\ Y \z /k, when Rf < R < R^. 
Using Eq. (18) and Eq. (19), V and f2 in the steady state are calculated as 



length of diffusion, / d = ^jD/a. Ay and \ z are the eigenvalues for the linearized equations 
for Y and Z, 




(24) 



(23) 




(25) 







(R < Rl R > R c n ) 




(26) 
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FIG. 2. (Color Online) Bifurcation of the motion of a droplet, (a) Phase diagram of the motion. 
One-point dashed line (blue) and dashed line (green) represent R = Id and R = Z a , respectively, (b) 
Schematic diagrams of the interface concentration field and the motion of a droplet in a laboratory 
frame. When the droplet moves straight (b-1) and rotates (b-2), the concentration field is mirror 
symmetric and asymmetric about the anterior-posterior axis, respectively, (c) Dependence of $7 on 
Lp 3 / (p 2 + l) 2 — 1. (d) Dependence of |0| on R. a = 1, = 1, and l a = 5. With these values, 
R\ = 0.83, = 4.55 and reaches a maximum value at R = 1.28. 



The interface concentration and the motion of the droplet in a laboratory frame are schemat- 
ically illustrated in Fig. 2 (b). When a droplet exhibits translational motion (R < R^,R> 
i?u), the interface concentration is symmetric about the anterior-posterior axis and the 
anterior-posterior axis always corresponds to the direction of the motion. On the other 
hand, when a droplet rotates (R\ < R < i?Jj), the interface concentration is asymmetric 
about the anterior-posterior axis and the particle is always the inside of the trajectory. Con- 
sidering the normalized Q, Cl = Q/Xy, Cl is the square root of the distance from the critical 
point, \Lp 3 / (p 2 + l) 2 — 1|, which corresponds to the stable solution of the normal form of 
pitch-fork bifurcation [10], as shown in Fig. 2 (c). Figure 2 (d) shows the dependence of \Q\ 
on R when la, l a and a are fixed. has a maximum value at a radius between R\ and 
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FIG. 3. (Color Online) Outline of the experimental system, (a) Schematic diagram. A petri dish 
with a radius of 7.5 cm was filled with 100 ml water. A droplet of oleic acid together with a 
3 mm particle of sodium oleate was placed on the water surface. The droplet moved spontaneously 
driven by the Marangoni effect [7]. (b) Trajectory of a droplet with a radius of 5.9 mm (300 (A). 
(b-1) Snapshots of a droplet per 1 s. The arrow indicates the particle fixed on the interface, (b-2) 
Trajectory of the center of mass of the droplet. Both scale bars represent 10 mm. 



V. COMPARISON WITH THE EXPERIMENT 



Using above theoretical results, we analyzed the corresponding experiment reported in [7] 
with additional data. The schematic diagram of the experimental setup is illustrated in 
Fig. 3 (a). One hundred milliliters of water, which was purified with a MilliQ filtering 
system (Millipore), was placed to a petri dish. A droplet of oleic acid (Wako Pure Chemical 
Industries; 159-00246) and a solid sodium oleate (soap) were floated on the aqueous phase 
in the petri dish. A solid column of sodium oleate (3 mm in length) was chosen from a 
commercial sample (Nacalai Tesque; 257-02). The movement of the solid/liquid composite 
was captured by a digital video camera at 30 frames per second at room temperature, and 
then analyzed using Image J (http://rsbweb.nih.gov/ij/docs/index.html). 

A time series of snapshots of a droplet, and the corresponding trajectory of the center 
of mass of the droplet are shown in Fig. 3 (b). Due to the fluctuation of rotational speed, 
the trajectory does not exhibit a closed circle. From the analysis of observed trajectories, 
the angular velocity, Q, was measured. Time courses of |0| are shown in Fig. 4 (a). Five 
experiments were performed for each volume to yield the distribution of which is shown 
in Fig. 4 (b). To determine the maximum value of the distribution, the distribution was fit- 
ted with the function, exp {— b A Q 2 (Q 2 — a) + c}, using the weighted nonlinear least-squares 
Levenberg-Marquardt algorithm, where the fitting parameters were a, b and c [11]. The 
number of events was taken as the weight for fitting. Using these parameters, the peak of 
the distribution, f2 typ , is calculated as when a < 0, and \Ja/2 when a > 0. We used 
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FIG. 4. (Color Online) Analysis of the experimental results, (a) Time series of the absolute value 
of angular velocity of a droplet with a radius of 2.3 mm (50 /A; a-1) and 5.6 mm (260 /A; a-2). 
(b) Distribution of |fi| obtained from experiments using a droplet with a radius of 5.6 mm (closed 
circles). Dashed line (orange) is the fitted curve. The parameters were estimated as a = 3.19, 
b = 0.78 and c = 6.48. One-point dashed line (blue) represents f2 typ = 1.3 rad/s. (c) Droplet 
size dependence of fityp- Gray circles represent O typ . Dashed line (orange) is the fitted curve. 
Estimated parameters were a = 5.4 s _1 , Id = 3.4 mm and Z a = 10.6 mm. 



Q t y P as the typical value of Q. To see the dependence of f2 tyP on the radius, f2 typ is plotted 
against the droplet radius in Fig. 4 (c). Here, the radius was estimated from the image of 
the droplet. As shown in Fig. 4 (c), the radius with the highest f2 typ and the critical radii, R\ 
and -R^, were all observed as predicted by our theoretical considerations. By fitting Eq. (26) 
to the obtained points, (R, fi typ ), we estimated 1^, l a and a. The estimated parameters 
of the fitted curve, which is shown in Fig. 4 (c), were a = 5.4 s _1 , Z a = 10.6 mm, and 
Id = 3.4 mm. It is noted that the time scale (the period of rotation) and the space scale 
(the size of the droplet and the particle) in the experiment are the same order as 1/a, and Id 
and Z a , respectively. So far, quantitative comparison is not possible within our model since 
a and (3 is difficult to measure experimentally. 
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VI. SUMMARY 



In summary, we manifested the physical mechanism of the spontaneous rotation of a 
mirror symmetric system consisting of a droplet and a particle attached at the interface. 
Using the flow field calculated with Stokes equation, we analytically found that when the 
Marangoni flow is strong enough, there are two critical radii of the droplet for rotation 
and the peak of angular velocity at a certain radius between these two critical radii. For 
verification, we compare our model with the experiment and these theoretical predicted 
results were ascertained in the experiment. Although we analyze only the spontaneous 
rotation driven by the Marangoni effect in this article, the advective transport of objects is 
generally observed in self-propelled phenomena such as droplets driven by any propulsive 
force and amoebae cells. Therefore, the advective nonlinearity similar to that in this article 
is expected to cause the symmetry breaking for rotation in various kinds of self-propelled 
phenomena. 
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